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Abstract 


This paper develops a real-time method based on the 
law of the wake for estimating forebody skin-friction 
coefficients. The incompressible law-of-the-wake 
equations are numerically integrated across the 
boundary layer depth to develop an engineering model 
that relates longitudinally averaged skin-friction 
coefficients to local boundary layer thickness. Solutions 
applicable to smooth surfaces with pressure gradients 
and rough surfaces with negligible pressure gradients are 
presented. Model accuracy is evaluated by comparing 
model predictions with previously measured flight data. 
This integral law procedure is beneficial in that skin- 
friction coefficients can be indirectly evaluated in real- 
time using a single boundary layer height measurement, 
hi this concept a reference pitot probe is inserted into the 
flow, well above the anticipated maximum thickness of 
file local boundary layer. Another probe is 
servomechanism-driven and floats within the boundary 
layer. A controller regulates the position of the floating 
probe. The measured servomechanism position of this 
second probe provides an indirect measurement of both 
local and longitudinally averaged skin friction. 
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Simulation results showing the performance of the 
control law for a noisy boundary layer are then 
presented. 
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A k control equation coefficient 
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longitudinal lateral in boundary layer, 
m/sec 

unsteady perturbation of lateral velocity, 
m/sec 

longitudinal coordinate, cm 

normal boundary layer coordinate, cm 

steady-state probe position, cm 

current probe position, cm 

commanded probe position, cm 

old probe position, cm 

normal displacement within laminar 
sublayer, cm 

probe position in boundary layer, cm 

point of transition between laminar and 
turbulent layers 

nondimensional normal boundary layer 

+ y f 

coordinate, y = --Re / -~ 

nonnal point of transition between laminar 
and turbulent layers 

lateral boundary layer coordinate, cm 

angle of attack, deg. 

Clauser nondimensional pressure gradient 
parameter 

intermittency factor 

differential pressure, kPa 

pressure error feedback, kPa 

frame interval, sec 

boundary layer thickness, cm 

boundary layer thickness, noisy 
simulation, cm 

estimate of boundary layer thickness, cm 

boundary layer displacement thickness, cm 

laminar contribution to boundary layer 
displacement thickness, cm 

turbulent contribution to boundary layer 
displacement thickness, cm 

controller damping ratio 

uniformly distributed random noise 
function 
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9 boundary layer momentum thickness, cm 

% 

9. m laminar contribution to boundary layer 

momentum thickness, cm 

* 

9 , b mrbulent contribution to boundary layer 

momentum thickness, cm 

K law of tire wake velocity distribution slope 

parameter 

surface roughness height, cm 

non-dimensional surface roughness, 

dummy variable of integration 

effective normal height of laminar 
sub-layer, cm 

p mean deviation in Cp 

If wake parameter 

<7 standard deviation in Cp 

oj () controller natural frequency, sec -1 



dp 

dx 

89 

8x 


longitudinal pressure gradient, kPa/m 
longitudhral momentum tlrickness gradient 


Mathematical Operators 


( ) first derivative with respect to tune 

( ) second derivative with respect to time 


Introduction 


Lifting body and wing body designs currently 
advocated for transatmospheric reusable launch vehicles 
(RLV) or for crew return from the space station are 
derived from variations on the original lifting body 
concept. 1 For a variety of reasons, these designs have 
large base areas relative to the vehicle size when 
compared to previous hypersonic designs. For example, 
on the X-33 and VentureStar™ (Lockheed Martin, 
Bethesda, Maryland) configurations the large base areas 
are required to accommodate aerospike engines. 
Conversely, the X-38 — derived from the original X-24A 
mold lines — has a large blunt base area resulting from 
the upper body flap required to trim the vehicle. In both 
of the above cases, the base area flow is highly separated 
resulting in large negative base pressure coefficients. 
Because of the large base-to-wetted-area ratios of these 
vehicles, the base drag makes up the majority of the 
overall vehicle drag. Thus the base drag has a large 
impact on critical vehicle performance parameters such 
as maximum payload and maximum cross-range. Any 


decrease in base drag has the potential to significantly 
improve the overall vehicle performance. An increase in 
the vehicle lift-to-drag ratio would allow for a much less 
severe glide slope during the terminal phase of the 
vehicle reentry. Reduction in glide-slope steepness 
would in turn make the terminal area energy 
management task considerably easier. 


Early Work Points to a Possible Solution 


There is a feasible means for decreasing base drag. 
For blunt-based objects with base areas that feature 
heavily separated flows, a clear relationship is 
demonstrated between the base drag and the viscous 
forebody drag. 2, 3 Figure 1 shows the base drag 
coefficient plotted against the viscous forebody drag 
coefficient. The viscous forebody drag coefficient is 
defined as the sum of all viscosity-induced drag 
components on the vehicle forebody. These components 
include skin friction, forebody separation, forebody 
wakes, and parasite drag due to protuberances. The 
trend presented in figure 1 shows that as the forebody 
drag is increased the base drag of the vehicle generally 
tends to decrease. This base-drag reduction is a result of 
boundary layer effects at the base of the vehicle. The 
shear layer caused by the free-stream flow rubbing 
against the dead, separated air in the base region acts as 
a “jet pump"’ 2 and serves to reduce the pressure 
coefficient in the base areas. The surface boundary layer 
acts as an insulator between the external flow and the 
dead air at the base. As the viscous forebody drag is 
increased, the boundary layer tlrickness at the aft end of 
the forebody increases — reducing the effectiveness of 
the pumping mechanism — and the base drag can 
potentially be reduced. 


Flight Tests on SR-71 LASRE Ext 


Recent work performed on the Linear Aerospike 
SR-7 1 Experiment (LASRE) flight program has 
demonstrated that the trend presented in figure 1 can be 
practically applied to reduce the base drag of a large 
projectile. 3 4 In the LASRE experiment, a portion of the 
forebody of the test model was roughened to increase 
the forebody drag. A decrease in base drag was observed 
at all Mach numbers — subsonic, transonic and 
supersonic. 


Trade-Off Between Forebodv Drag and Base Drag 


More significantly, figure 1 data show that projectiles 
with a base drag coefficient greater than 9.39 
(referenced to the base area) must necessarily lie on the 
relatively steep segment of the curve presented. In this 
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region a small change in the forebody friction drag 
should result in a relatively large change in the base 
drag. Conceptually, if the added change in forebody skin 
drag is optimized with respect to the base drag, then it 
may be possible to reduce the overall drag of the 
configuration. That is, based on Hoerner’s data, an 
optimal drag bucket should exist. Figure 2 shows this 
optimal drag region. Here the base drag and the total 
vehicle drag are plotted against the viscous forebody 
drag including measured data for several hypersonic 
lifting-body and wing-body configurations: 3 These are 
the X-15, M2-F1, M2-F2, Shuttle, HL-10, X-24A, 
X-24B, and the SR-71 LASRE. The simple model in 
figure 2 is presented only as an illustration of the drag- 
bucket concept. Clearly, the forebody pressure profile 
and the presence of induced drag and localized 
interference or flow separation drag will likely alter the 
shape of the optimal curve presented in figure 2. The 
challenge is to modify the forebody drag coefficient to 
produce the minimum overall vehicle drag in a real- 
world configuration. 

On the Need for Real-Time Forebody Viscous Drag 
Measurements 

An intriguing concept for implementing the drag 
optimization described above utilizes recent advances in 
microelectromechanical systems (MEMS) 4 5 technology. 
In this concept MEMS technology is used to adaptively 
vary the surface roughness in order to track the optimal 
skin friction. This design uses a “fish-scale” pattern of 
microscale actuators distributed along the vehicle skin. 
The position of these microscale surface panels is 
commanded in real-time to modify the local surface 
roughness and increase or decrease the total forebody 
viscous drag coefficient, as required, to achieve the 
overall drag minimum. This optimization process 
requires that the viscous forebody drag coefficient be 
sensed and available for feedback to the MEMS 
controller in real-time. To be practical for the 
optimization problem, the forebody drag measurement 
device must be self-contained; that is, the measurement 
system must be installed entirely onboard the vehicle. 
Additionally, the device must be able to adapt to widely 
varying boundary layer thickness. 

None of the existing methods for sensing the surface 
skin friction are easily adaptable to the above 
optimization problem. This paper briefly examines why 
there are limitations with the current skin-friction 
measurement techniques and then proposes a new 
method for sensing the viscous forebody drag 
coefficient based on the law of the wake. 6 The 
engineering model developed for this new method 


allows longitudinally averaged skin-friction 
measurements to be made in a direct maimer with a 
single measurement of the boundary layer thickness 
near the aft end of the forebody. For simplicity, all 
analyses and data presented assume incompressible flow 
relationships. It is simply a matter of detail to extend 
these discussions to account for the effects of 
compressible flow. 

Note that use of trade names or names of 
manufacturers in this document does not constitute an 
official endorsement of such products or manufacturers, 
either expressed or implied, by the National Aeronautics 
and Space Administration. 

Background 

This section briefly discusses existing methods for 
sensing surface skin friction and presents the reasons 
why these methods are unsatisfactory for our current 
measurement problem. Current methods for measuring 
surface skin friction can be grouped into two classes: 
direct measurements and indirect measurements. 

Direct Skin Friction Measurements 

The direct measurement methods sense surface 
shearing stress by using small floating balances or 
through the use of surface-imaging interferometry. The 
balance methods are generally self-contained and have 
been successfully applied to flight vehicles. 7,8 
Unfortunately, experience has demonstrated that 
surface-mount skin-friction gauges are also susceptible 
to damage and thermal shifts, and they require frequent 
re-calibration to achieve accurate results. 

Oil-film interferometry has been used to produce vety 
accurate localized skin-friction measurements; 9 
however, these techniques are usually applied only in a 
laboratory setting. Because they require remote optical 
imaging of the flow field, interferometry methods are 
inherently not self-contained and are not easily adapted 
for use on flying vehicles. The rapid development of 
fiber-optic sensors may provide a practical method to 
overcome this problem. Developing these optically 
sensed skin-shear sensors to be lugged and reliable will 
be a challenge for some time yet to come. 

Indirect Skin Friction Measurements 

Two of the most common indirect methods used to 
compute skin-friction coefficients on flight vehicles are 
Preston tubes and boundaiy -layer rakes. Both types of 
sensors are reasonably easy to fabricate and install, and 
offer the advantage of being self-contained, reliable, and 
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rugged. Unfortunately, both designs have limitations 
that restrict their real-tune applicability and make them 
generally unsuitable for the drag optimization problem. 

Preston Tubes 10, 11 offer a simple hardware solution 
where a single pitot tube is located as close to the 
surface as possible. By assuming a log-law velocity 
profile, 12 and using measurements of the local wall 
static and surface pitot pressure, the local skin-friction 
coefficient can be calculated. These methods have been 
successfully applied to a large number of flight and 
wind tunnel tests. However, the drag optimization 
problem requires that the entire viscous forebody drag 
coefficient — which includes the effects of skin friction, 
forebody separation, forebody wakes, and parasite drag 
due to protuberances — must be sensed. Unfortunately, 
the log-law velocity distribution is not valid in the outer 
layers of the boundary layer or within the wake region. 
Thus it is believed that Preston tubes will not be able to 
completely capture the range of required forebody drag 
measurements. It is possible that tire Preston tube 
analyses can be modified to allow for tire effects of the 
outer wake region; however that analysis is beyond the 
scope of this paper. 


A New Measurement Approach 

The real-time measurement technique being 
developed in this paper exploits design features of both 
Preston tubes and boundary-layer rakes. In this 
approach a reference pitot probe is inserted into the 
flow, well above the maximum anticipated thickness of 
the local boundary layer. Another pitot probe is mounted 
on a servomechanism which moves in a direction 
normal to the surface. A controller regulates the position 
of the servomechanism based on the measured pressure 
difference between the moving pitot probe and the 
reference pitot probe. This design concept is pictured 
schematically hi figure 3. The design will be tailored to 
typically allow the probe to track the area within the 
outer wake region near the edge of the boundary layer. 
A position sensor on the probe gives a measure of the 
local boundary layer thickness. This probe position is 
related to the forebody skin-friction coefficient using a 
law-of-wake analysis. This analysis is presented in the 
next section. The ciux of this analysis is the 
development of an integral relationship that describes 
the skin-friction coefficient in terms of a local 
measurement of the boundary layer thickness. 


Boundary- layer rakes consist of an assembly with 
multiple pitot probes placed in a normal line away from 
the skin of the vehicle. Boundary-layer rakes sense the 
local dynamic pressure profile within the boundary layer 
and use these data to compute the velocity distribution. 
This velocity distribution is numerically integrated to 
produce a measurement of the local momentum 
thickness. 12 


Integral Law of the Wake Analysis 

The key to the floating-probe concept described in the 
previous section is Coles’ law of the wake. 6 The law of 
the wake relates the outer (turbulent) boundary-layer 
variables to the nondimensionalized boundary-layer 
coordinates 


In order to obtain accurate momentum thickness 
measurements using a boundary-layer rake, it is critical 
that the rake be designed to obtain adequate normal 
resolution within the critical shearing layers near the 
surface where the momentum defect distribution 
reaches a peak value. This “hump” hi the momentum 
thickness curve occurs near the interface of the 
laminar sublayer with the outer turbulent log-law 
velocity distribution layer. This resolution requirement 
means that a large number of measurements, typically 
6-10 probes, must be made in order to sense the 
velocity profile with a sufficient amount of resolution. 
This large amount of data must be sampled and 
integrated in real time. Additionally, the frame-by-frame 
results must be heavily post-processed to eliminate the 
effects of local boundary noise. While these limitations 
certainly do not preclude the use of boundary-layer 
rakes for real-tune forebody drag measurements, it is 
believed that these results can be achieved by a more 
straightforward measurement technique. 


+ B (1) 

where, 



8 is the local boundary-layer thickness, K is the 
slope parameter, B is a bias parameter, n is the wake 
parameter, and C F is the local viscous forebody drag 
coefficient. The currently accepted “best value” for K is 
approximately 0.41. The bias parameter, B, varies with 
the level of surface roughness and for a smooth plate has 
a numerical value of approximately 5.0. The wake 
parameter is proportional to the local longitudinal 
pressure gradient, and for a zero-gradient flow condition 
this wake parameter has a numerical value of 
approximately 0.55. The roughness-dependent bias term 


+ i 

u 
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ln[y + ] + 2nsin 2 
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can be eliminated from equation 1 by expressing the law 
of the wake in terms of the local velocity defect. 


1 - 


u(yY 

u. . 



2IIcos2 


“7t y~ 
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General Solution Method for Local Skin-Friction 
Coefficient 

As developed in Appendix A, equation 3 is integrated 
across the depth of the boundary layer, including the 
laminar sublayer, to relate the momentum and 
displacement thickness to the local boundary thickness 
and the local pressure gradient. The resulting 
relationships are (1) momentum thickness and 
(2) displacement thickness. 

(1) Momentum thickness: 




[0.16425 + 2.66116V0.88881 
[[10.4874 + 10.6227 p + 2.13276^0.88881 + p] Jcj 


55.9194- 40.6563 VO.88881 + p 


Re 


(4(a)) 


Jcj l- 137.422 + 162.4 p + 32.605 W0.88881 + p] 


Re 


gradients with ||3| >0.1 is considered to be strong. 
White 14 develops a similar set of integral relationships 
for the momentum and displacement thickness, but does 
not include the effect of the laminar sublayer. 

Closure is added to the problem by applying the 
correlation equation developed by White, 14 


c 


f ~ 


o.3e 


lo g io Re x ? 


1.338*/e 


-,1.74 + 0.31 8Ve 


(6) 


Equations 4, 5, and 6 can be used to solve for the 
entire state of the local boundary in terms of the 
normalized boundary-layer thickness, Reynolds 
number, and the Clauser pressure gradient parameter. 
Parameters which are estimated include the local skin- 
friction coefficient, momentum thickness, displacement 

thickness, and the boundaiy-layer shape parameter 

* 

(8 /0). This system of equations is extremely 
nonlinear but is reasonably well-behaved and is 
numerically solved using a simple substitution iteration 
procedure. Unfortunately, the correlation of equation 6 
was developed from smooth surface data only. Thus this 
solution procedure should be applied to smooth surfaces 
only. For the remainder of this paper, this iterative 
solution procedure will be referred to as the nonlinear 
method. 


(2) Displacement thickness: 

* 

— = ~Jc7[[ 0.16425 + 2.66116^0.88881 + P]] 

X X (4(b)) 

55.9194 - 40,6563 VO, 88881 + (3 


Zero Pressure Gradient Solution for Longitudinally 
Averaged Skin-Friction Coefficient 

If the local pressure gradient is small or negligible 
(P~0), then the longitudinally averaged forebody 
viscous drag coefficient is directly related to the local 
momentum thickness using the integral fonn of the 
Von Karman momentum equation 


hi equations 4(a) and 4(b), P is the Clauser 
parameter, and is related to the local pressure gradient 
by 


P = — — (5) 

C f-q 

Clearly negative values of P correspond to favorable 
pressure gradients and positive values of P correspond 
to adverse pressure gradients. A pressure gradient with 
IPI <0.1 is considered to be a weak, and a pressure 


1 pr c _f 
xJo 2 


dx 



(7) 


Substituting equation 7 into equation 4(a) and 
evaluating for P = 0 , an expression for the 
longitudinally averaged viscous drag coefficient in 
terms of the local skin-friction coefficient is derived 


C F = ^JT f [5.3464-24.9987^] 
+ -^-[35. 12- 213.342 TcT] 


(8(a)) 
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Ill equation 8(a), C p is the longitudinally averaged 
viscous drag coefficient. Evaluating equation 1 at 
y = 8 gives an implicit formula for the local skin- 
friction coefficient in terms of the 8 and Re x 


— = 2.439 


In 


C f 5 i? 

2 x Re * 


+ 2n 


+ B (8(b)) 


As mentioned earlier, for a smooth surface (1 = 5.0 . 
For rough surfaces, White 14 has curve-fit the 
Prandtl-Schlicting sand grain roughness cur ve to obtain 
the following formula 


B = 5.0-2.4391n[l +0.30 k + ] (9) 

In equation 9, K ' is the nondimensional surface 
roughness and is defined as 



A surface is defined as “rough” when k' > 5.0 . For 
k < 5.0 , the surface is said to be “hydraulically 
smooth.” Equation 9 is valid for turbulent flow 
conditions that are smooth, transitionally rough (only 
part of the surface is rough), and fully rough (the entire 
surface is rough). Collectively, equations 8, 9, and 10 
offer a nonlinear set of equations which can be used to 
solve for the longitudinally averaged viscous drag 
coefficient in terms of the local boundary-layer 
thickness, Reynolds number, and the surface roughness 
ratio. This set of equations is well-behaved and can be 
solved iteratively using a simple substitution method. 
Since the local skin-friction coefficient is computed by 
evaluating the law of the wake at the edge of the 
boundary layer for the remainder of this paper, this 
solution method is referred to as the edge method. 


Closed-Form Solution Method for Averaged Skin- 
Friction Coefficient 

When the entire surface is rough (a fully rough flow 
regime), experimental observations show that the ratio 
of file local and averaged skin-friction coefficients is a 
function of only the normalized surface roughness, 
(x/k s ). For a fully rough flat plate Mills, 15 et. al. have 
developed the following relationship for the skin- 
friction ratio 



2.635 + 0.618 In 



3.476 + 0.707 In 



(ID 


Equation 1 1 is valid over a Reynolds number range 

c o 

from approximately 10 to 10 . For large Reynolds 

n 

numbers (Re x > 10 ) an approximate skin-friction ratio 
for fully developed turbulent flow on a hydraulically 
smooth plate is derived by setting the roughness ratio to 
a very large value. Evaluating equation 11 at 
x/K s = 10 8 , gives G[x/k s ] ~ 0.896. 

By substituting G[x] into equation 8(a) and 
simplifying, the solution for the averaged viscous drag 
coefficient are written in closed form 

Equation 12 gives the relationship between the local 
boundary layer thickness and the averaged surface skin- 
friction coefficient. This simple result is important 
because it requires only a single measurement of the 
boundary layer thickness and an estimate of the 
Reynolds number and surface-roughness ratio to 
determine the local momentum defect. When the 
pressure gradients are small, the local momentum defect 
is a good approximation of the averaged skin friction. 
Data to be presented in the Results and Discussion 
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section demonstrate that the closed-form model gives 
results that ate accurate to within 5 percent as long as 
IPI < 0.1 . This accuracy estimate holds for fully rough 
surfaces with turbulent flow and for hydraulically 
smooth surfaces with turbulent flow occurring at the 
leading edge. For the remainder of this paper 
equation 12 is referred to as the closed-form solution. 

Calculations from the various solutions methods are 
compared in the Results and Discussion section. These 
comparisons help in determining the regions of 
applicability for the solution methods. 

Design of the Boundary Layer Tracking Probe 



(14(a)) 


(14(b)) 


As mentioned earlier the integral form of the law of 
the wake — where the skin-friction coefficient is 
expressed in terms of the local boundary layer thickness 
— suggests an innovative concept for a skin-friction 
measurement device. While the actual mechanical 
design of the probe can vary substantially, the basic 
concept is that a reference pitot probe is inserted into the 
flow, well above the maximum anticipated thickness of 
the local boundary layer. Another pitot probe is 
servomechanism-driven to move normally to the 
surface. A controller regulates the position of the 
floating probe. As developed in Appendix B the 
continuous-time control law for this system is 

y + 2 C°v + % 2 y = “» 2ifr q (13) 


In equation 14 k is the time-frame index and At is the 
time interval between each frame. These equations can 
be implemented recursively starting with the initial 
position of the tracking probe. Since the measured 
dynamic pressure difference is a function of the local 
position within the boundary layer, this control law is 
clearly nonlinear. As such, it is expected that the control 
law behavior will differ somewhat from the response of 
a linear second-order system. Clearly though, n is the 
control gam parameter. An example probe control-law 
response for an unsteady boundary layer will be 
presented in the Results and Discussion section. 

Final Probe Position 

When the probe is positioned within the boundary 
layer so that 


In equation 13, y is the probe position measured 
relative the to the surface, q is the local dynamic 
pressure at position y within the boundary layer, and 
Aq is the difference between the dynamic pressure at 
the edge of the boundary layer and the local dynamic 
pressure. Choosing the natural frequency (ot>„) and 
damping ratio (^) so that the open- loop system 
eigenvalues are stable, causes the unforced probe 
position to always decay to zero. The decaying tendency 
of the open-loop system insures that the measurement- 
probe will not remain stuck outside of the boundary 
layer where the pressure difference between the 
reference and measurement-probes is zero. The analog 
flow diagram for this control equation is depicted in 
figure 4. 

Discrete Control Law 

Equation 13 is discretized using the trapezoidal 16 
rale. After performing some extensive algebra, the 
discrete control equation is written as 


2L q J k 


= 0 


equation 14 is rewritten as 


(15) 


yt+i-yk 


1 + 2C%f 




t] 


(16) 


Since C , , co /; , and 8/ are all positive constants, 
equation 16 implies that 

K+t-^I^K-^-tl < 17 ) 

Thus file difference equations converge to a steady 
value. As a result the probe seeks a final position where 


'Acyl _ 2 

- q - final w 


( 18 ) 
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Equation 18 is written in terms of the final velocity 
ratio as 


' u “ 

_ 

r<?i 

~ 1 ” 

u 

e - 1 

final W 

lq e \ 

final ^ n + 2 


(19) 


Equation 19 is used to correct the final probe position 
by giving the proper estimate of the boundary layer 
height, assuming that the velocity profile in the outer 
layers of the boundary layer are approximated by a 
power-law function. The resulting correction factor is 


y final 


f n + 21"'/ 2 


( 20 ) 


In equation 20 the parameter rn is the velocity profile 
exponent and is the function of the local Reynolds 
number. Values for m, which are required to perform the 
correction of equation 20, are presented in the Results 
and Discussion section. Once the final probe position 
has been achieved; the boundary layer thickness 
estimated by equation 20 is used to evaluate the local 
and averaged skin-friction coefficients using whichever 
of the three solutions methods (developed earlier) is 
appropriate. 


Results and Discussion 

This section is divided into two parts: a Model 
Validation section and a Probe Response section. In the 
Model Validation section the three solution methods for 
the law of the wake model are compared. These 
comparisons establish the errors that pressure gradients 
cause in the edge and closed-form solution methods. 
Next the closed-fonn solution is compared to the edge 
solution for a variety of surface roughness. Finally, skin 
friction coefficients and velocity profiles predicted by 
file law-of-the-wake model are compared to flight- 
obtained boundary layer measurements. In the Probe 
Response section an example of dynamic response of 
the probe control law for a noisy boundary layer is 
presented. 

Model Validation 

First an estimate of errors in the edge and closed-form 
solution methods caused by longitudinal pressure 
gradients are examined. These results are presented in 
figure 5. In figure 5 the top graph shows the local skin- 
friction coefficient predicted by the nonlinear model as a 
function of Reynolds number. The curves correspond to 
values of p = {-0.25,0.0,0.25, and 1.0}. The 
solutions for the edge and closed-form methods ( P = 0 


implicit) are also plotted. For the closed-form solution 
the “smooth-plate” value of G[x/k s ] = 0.896 was 
used. For P = 0 the nonlinear and edge solution 
methods give nearly identical results. The closed-form 
solution shows approximately a 2-percent error at 
Re x = 10 . This error diminishes to nearly zero at 
Re x = 10 8 . The bottom graph in figure 5 shows the 
percent deviation in the closed-form model, defined by 

c f -c f 

percent of deviation = 100 x — (21) 

c f 

J nonlinear 

This parameter is plotted as a function of Reynolds 
number for values of p = {-0.25, 0.0, 0.25, and 1.0}. 
The calculated deviations in the edge method are 
smaller than those shown in the figure 5 bottom graph 
by approximately 1 percent. These data are not shown. 
The deviations shown in tins bottom graph of figure 5 
are also representative of the errors that are expected for 
the averaged skin friction values. Clearly, there is error 
introduced by the simplifying assumptions used in 
deriving the closed-form model. However, if the local 
pressure gradient is small (|p| <0.10), the errors are 
less titan 5 percent over a very wide range of Reynolds 
numbers. 

The simplifying assumption was introduced into the 
closed-form model by the simplifying assumption of 
equation 1 1 . However, as shown in figure 5, this error is 
less than 2 percent over a wide Reynolds number range 
for a hydraulically smooth surface. Figure 6 shows how 
this error varies as a function of surface roughness. The 
figure 6 top graph compares the edge solution for 
averaged skin-friction coefficient to the closed-form 
averaged skin-friction coefficient. The data presented on 
this curve were computed for a Reynolds number of 
5.0 x 10 with roughness levels varying from a very 
smooth surface to a vety rough surface. Also plotted on 
this figure 6 top graph is the transitionally rough skin 
friction roughness formula derived by White. The 
level of agreement is very close. The bottom graph of 
figure 6 shows the percent deviation plots (where the 
edge solution is assumed to be the “truth set”). The 
deviations are all within 2.5 percent, even at the low 
Reynolds numbers. As the Reynolds number and surface 
roughness increase, the deviations diminish to 
approximately zero. 

Comparisons to Flight Data 

Figures 5 and 6 show data in the absence of 
significant longitudinal surface pressure gradients. In 
these two figures the closed-form model predicts the 
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averaged skill friction to within 5 percent when 
compared to the more complex iterative solution 
methods. This result justifies using the simple model of 
equation 12 for the real-time measurement system. 
Further justification is evident in comparing averaged 
skin friction predictions of the closed-form model 
against data obtained from two flight experiments 
conducted by Saltzman and Fisher. 17, 18 During these 
series of tests, both supersonic and subsonic data were 
obtained; however, since this analysis is currently 
applied only for incompressible flow conditions, only 
the subsonic flight data were analyzed. In these flight 
experiments, boundary-layer rake measurements were 
obtained on the underside to the fuselage of an A-5A 
aircraft, and on the whig of the XB-70-1 aircraft. For 
both series of flight tests a calibrated Preston tube was 
also installed near the rake to supply additional 
measurements of the local skin-friction coefficients. 
Figure 7 shows the A-5A instrument locations and 
figure 8 shows the XB-70-1 rake and Preston tube 
locations. Table 1 summarizes the data points that were 
analyzed. 


Table 1. Summary of subsonic flight data analyzed. 


Data 

point 

Vehicle 

Re x 

a, deg 

^ edge 

x , cm 

1 

A5A 

1.54e+7 

6.2 

0.82 

508.0 

2 

A5A 

2.20e+7 

3.8 

0.64 

508.0 

3 

A5A 

3.00e+7 

6.4 

0.60 

508.0 

4 

A5A 

3.1 le+7 

7.0 

0.51 

508.0 

5 

A5A 

4.94e+7 

3.2 

0.81 

508.0 

6 

A5A 

7.42e+7 

0.3 

0.81 

736.0 

7 

XB-70-1 

1.14e+8 

8.7 

0.41 

1557.5 

8 

XB-70-1 

1.18e+8 

8.1 

0.49 

1557.5 

9 

XB-70-1 

1.24e+8 

8.2 

0.46 

1557.5 

10 

XB-70-1 

1.30e+8 

8.9 

0.53 

1557.5 

11 

XB-70-1 

1.33e+8 

8.3 

0.51 

1557.5 

12 

XB-70-1 

1.41e+8 

5.5 

0.94 

1557.5 

13 

XB-70-1 

1.57e+8 

5.2 

0.95 

1557.5 

14 

XB-70-1 

1.73e+8 

5.7 

0.91 

1557.5 


For the A-5A tests, all of the data points analyzed 
were at angles of attack greater than 3.8°, except point 
number 6 for which the angle of attack was 0.3°. 
Because the boundary-layer rake for these tests was 
located on the bottom of the vehicle, for the higher 
angles of attack (points 1—5) the effect of the noseboom 


on the Reynolds number was ignored. For these points a 
total run length of 508 cm was used. For the low angle- 
of-attack data point (point 6) the effect of the noseboom 
became significant and a run length of 736 cm was used. 
For the XB-70-1 tests the run length was approximately 
1557.5 cm. 

For each data point the boundary layer thickness was 
obtained by fitting the data with a 1/m' 1 ' power law 


lny 


hi 8 + min 


jj-/ 0.99 


( 22 ) 


Figure 9 shows an example linear-log plot for data 

point number 2 (A-5A) in table 1 . The y-intercept value 

is 1.8112 and the boundary layer thickness is 
18112 J * 

approximately 8~e ' =6.118 cm. The linear- log 

slope represents the power law exponent, hence 
m ~ 7.9613 . Based on this “measured” boundary layer 
thickness, the closed-form solution was used to compute 
both the averaged and local skin-friction coefficients. 
These data were then used with equation 3 to compute 
the velocity profile and momentum distribution within 
the boundary layer. Figure 10 shows the velocity and 
momentum profile comparisons for data point 2 (A-5A, 
Re x = 2.20 x 1 0 7 ) . Outside of the laminar sublayer, the 
comparisons are excellent. Similar velocity and 
momentum distribution profiles are presented in 
figure 1 1 for data point 10 (XB-70-1, Re x = 13.0 x 10 7 ). 
The comparisons are very similar to those in figure 10. 
The characteristics of these selected data points are very 
similar to all of the data points which were examined. 

Figure 12 summarizes the flight data analysis. 
Figure 12 plots the law of the wake (averaged) skin- 
friction (closed-form solution) coefficients calculated 
using flight-measured boundary layer thickness against 
previously published ’ values detenu ined by the 
authors. The published data include local skin friction 
estimates deteimined from the rake data using an 
adaptation of the Clauser method. Flight -measured 

Preston tube data are also shown. The Preston tube data 
were analyzed using the log-law calibration method of 
Hopkins and Keener. 10 The published data were 
expressed in teims of the local skin-friction coefficients 
and must be transfoimed to averaged values for the 
comparisons to be valid. The local skin friction data 
were transfoimed to averaged values using the method 
of Mills. 1 '’ In this method the skin friction ratio is 
derived from integrating the local skin friction from the 
leading edge of the forebody along a flow streamline to 
the local position. The effect of laminar flow near the 
leading edge is accounted for. The result is 
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C F 

C f 


2.9187, , 

r^j 

In [0 06Re ] 

J R %, 

L Re xJ 


+ 1.1495 


Re x hi 2 [0.06Re x ] ' 
Re x In 2 [0.06Re ] 

tr X (r J 


calculations directly by simply measuring the local 
boundary layer thickness. This section analyzes the 
dynamics of the tracking device that is used to obtain 
(23) the boundary layer thickness measurement. In this 
section the performance of the tracking probe is 
modeled using both steady-state and noisy boundary 
layer simulations. 


In equation 23 Re is the transition Reynolds 

ir 

number and Re x is the local Reynolds number. For this 
transformation a transition Reynolds number of 
approximately 5.0 X 10 5 was assumed. 

Figure 12 shows the skin-friction coefficients plotted 
against the normalized boundary layer height. The 
(adjusted) published data generally scatter about the 
law-of-tlie-wake curve; however the scatter is large. 
Table 2 summarizes the mean and standard deviations of 
the adjusted rake and Preston tube Cp values away from 
the law of the wake determined values. 


Table 2. Deviations of flight measured 
C F from law of the wake. 



Rake 

Preston tube 


F, 

percent 

O', 

percent 

F, 

percent 

o, 

percent 

A-5A 

10.68 

31.0 

1.2 

27.0 

XB-70-1 

-5.9 

11.7 

-5.3 

10.7 


The XB-70-1 data show a much better agreement than 
do the A-5A data. The reasons for the large amount of 
scatter in the A-5A data ate unclear. Because the test 
array was located on the underside of the vehicle 
fuselage, it is possible that angle-of-attack effects 
induced strong local pressure gradients which had a 
significant influence on the local skin friction 
measurements. The 5- to 10-percent mean deviations are 
not, however, unreasonable. 

Simulation of the Boundarv-Laver Tracking Probe 
Response 

The previous section established that the integral law- 
of-the-wake analysis is a reasonably accurate predictor 
(within 5 to 10 percent) of the average skin-friction 
coefficient in the absence of large external pressure 
gradients. For real-tune applications the closed-fonn 
solution given by equation 12 is the preferred solution 
method. The value of the closed-fonn solution is that the 
averaged skin-friction coefficient can be incorporated in 


Effects of the Control Gam Parameter 

In this series of simulations the boundary layer 
thickness and the velocity profile are assumed to be 
constant as a function of time. The probe is released at 
time t = 0 with a starting value of y/8 = 0.25 . To 
match the expected transonic Reynolds numbers that 
atmospheric reusable launch vehicles encounter, the 
Reynolds number was held constant at 

7 

Re x = 5.0 X 10 . The simulations were all run at 
250 samples per second. Figure 13 demonstrates the 
effect of file control gain parameter (n) on the system 
response. Here the simulation is run with the natural 
frequency, damping ratio, and Reynolds number fixed 
at co n = 4k sec \ C, = 0.7071, aud Re x = 5.0 X 10 ? 
The control gain parameter is varied from 10 to 500. 
The results illustrate an interesting point regarding the 
trade-off between final position of the probe and the 
oscillatory stiffness of the model response. As predicted 
by equation 19, when the gain is increased, the final 
position of the probe is closer to the boundary layer 
edge. Unfortunately, high gain values induce an 
undesirable response in the system, and the probe never 
settles down to a stable value. Clearly, it is preferable to 
allow the probe to track at some fraction of the 
boundary layer height, and correct the final position 
using equation 20. It is recommended that values of n be 
kept somewhere within the range from 20 to 50. 

Correction for the Final Probe Position 

The velocity profile exponent (m), which was required 
to correct the final probe position and estimate file true 
boundary layer thickness, was determined by running a 
series of simulations in which the Reynolds number aud 
gain parameter were allowed to vary. These results are 
presented in figure 14. Figure 14 shows the variation of 
m with Reynolds number for control gains of 10, 20, 50, 
100, and 200. The data of figure 14 were fit with a 
simple linear model of the form 

m[n, Re x ] = C Q [«] + C, [»]log ^[ReJ (24) 
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The coefficients of equation 24 are displayed in 
Table 3. 


Table 3. Velocity profile exponent: fit coefficients. 


n 

C 0 

Ci 

10 

-0.51785 

0.95893 

20 

0.11626 

0.96786 

50 

0.69694 

1.06950 

100 

1.29530 

1.14760 

200 

1.89040 

1.24180 


Linear interpolations of the coefficients in Table 2, as 
a function of n, are used to obtain a correction function 
for the final probe position. 

Tracking an Unsteady Boundary Laver 

The steady-state analyses of the previous section were 
extremely useful in understanding the dynamics of the 
tracking probe control law; however, a true boundary 
layer is not at all stationary, hi reality a turbulent 
boundary layer has random elements and instantaneous 
variations that can only be statistically predicted. The 
result is an intermittent boundary layer with a ragged 
shape. In the outer regions of the boundary layer there 
are turbulent fluctuations which extend beyond the 
normally defined boundary layer edge. Conversely, 
there are intermittent regions inside of the normally 
defined boundary layer where the velocity is equal to the 
free-stream velocity. The law-of-the-wake analysis 
presented in this paper deals only with the time- 
averaged velocity components of the boundary layer. 
However, the dynamic flow components will effect the 
performance of the boundary layer tracking probe and 
this flow intermittency must be modeled in the probe 
simulation. In this paper the random components within 
the boundary layer are modeled using data taken from 
Klebanoff. 19 This intermittent boundary-layer model is 
developed in detail in Appendix B. 

The response of the probe to an intermittent boundary 
layer will now be presented. Reynolds number was set 
at Re x = 5.0 x 10 7 , and the natural frequency, damping 
ratio, and control gain parameter were varied to 
determine the best combinations required in order to 
accurately track the edge of the boundary layer. Using 
the corrected probe position, the skin-friction coefficient 
is computed using equation 12. Figure 15 presents the 
results of a noisy simulation run with a very stiff control 
law. In this simulation the parameters were set at 


co n = ji/2sec _1 , £ = 0.7071, and n = 25. The 
response of the tracking probe and the corrected probe 
position are shown relative to the normalized boundary 
layer edge hi figure 15(a). The estimated skin-friction 
coefficient calculated using equation 12 is shown 
compared to the known value in figure 15(b). The probe 
position correction based on equation 25 works quite 
well and the estimated skin-friction coefficient tracks 
within 5 percent of the “truth value”-once the initial 
startup transient has died out. This simulation is 
presented for illustrative purposes only. The particular 
values which are preferred for the control law 
parameters must be tuned to the test conditions, 
hardware characteristics, desired probe response time, 
and noise rejection characteristics. Clearly more work 
needs to be performed in this area. 

Summary and Concluding Remarks 

Lifting body and whig body designs currently 
advocated for transatmospheric reusable launch vehicles 
are derived from variations on the original lifting body 
concept. These designs have large base areas relative to 
tlie vehicle wetted area. Thus the base drag has a large 
impact on critical vehicle performance parameters such 
as maximum payload and maximum cross-range. Any 
decrease in base drag can have the potential to 
significantly improve the overall vehicle performance. 
Recent work has demonstrated that by increasing the 
forebody drag of these vehicle shapes it may be feasible 
to reduce the base drag so significantly, that the overall 
drag of the vehicle can be reduced. Methods for 
optimizing this drag reduction process are currently 
being studied. 

Any attempt to implement the drag optimization 
requires a practical method for measuring the viscous 
forebody drag of the vehicle in real time. To be practical 
for the optimization problem, skin friction 
measurements must be self-contained, that is, they must 
be installed onboar d the vehicle entirely. As discussed, 
none of the existing methods for sensing the surface 
skin friction are readily adaptable to the above 
optimization problem. 

This paper proposes a new method based on the law 
of the wake. In this concept a reference pitot probe is 
inserted into the flow, well above the anticipated 
maximum thickness of the local boundary layer. 
Another probe is servomechanism-driven and floats 
within the boundary layer. A controller regulates 
position of the floating probe. The measured 
servomechanism-position gives a direct measurement of 
tlie local boundary layer height. The measured boundary 
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layer height is related to the local and longitudinally 
averaged skin friction by an integral form of the law of 
the wake. 

In the law-of-tlie-wake analysis the turbulent 
boundary layer velocity profile equations were 
numerically integrated across the boundary layer depth 
to develop an engineering model that relates 
longitudinally averaged skin-friction coefficients to 
local boundary layer thickness. The law-of-the-wake 
analysis develops three solution methods: 1) an iterative 
nonlinear solution, 2) an edge solution where the local 
skin friction is evaluated from the law of the wake at the 
edge of the boundary layer, and 3) a less-general closed- 
form solution. The non-linear solution method solves 
for the local skin friction coefficient, momentum 
thickness, displacement thickness, and the boundary 
layer shape parameter. The nonlinear solution method is 
generally applicable to smooth surfaces with arbitrary 
longitudinal pressure gradients, but does not allow the 
averaged skin-friction coefficient to be evaluated. The 
edge method is valid only when the pressure gradients 
are small, but allows both the local and averaged skin- 
friction coefficients to be solved for. The edge method is 
valid for smooth, transitionally rough, and fully rough 
surfaces. The closed-form solution method is the least 
general of the tluee solution methods and is applicable 
only for fully rough surfaces when tire pressure gradient 
is small. The closed-fonn solution offers the cleat- 
advantage of being an explicit solution; whereas, the 
first two methods require an iterative calculation. For 
real-time applications, as long as the flow constraints are 
not violated, the closed-form method is clearly the 
preferable choice. 

Comparisons of the closed-form solution method to 
flight data show that the velocity longitudinally 


averaged skin friction can be sensed with reasonable 
accuracy. The flight data show that it is clearly possible 
to achieve mean accuracies on the order of 5 to 
10 percent. Applying corrections for the longitudinal 
pressure gradient and the effects of compressibility can 
likely reduce these errors further. 

The dynamics of the tracking device used to obtain 
the boundary layer thickness measurement were 
presented. The performance of the tracking probe was 
modeled using both steady-state and noisy boundary 
layer simulations. Simulation results showing 
performance of the boundary layer tracking probe for a 
noisy boundary layer were presented. The particular 
values that are preferred for the control law parameters 
must be tuned to the test conditions to obtain the desired 
probe response tune and noise rejection by the probe 
control law. 

To be more generally applicable, additional 
comparisons of the integral law of the wake to transonic 
and supersonic flight data should be performed. If 
possible, a compressibility correction must be 
developed. Additional data comparisons should also be 
performed for known pressure gradients, a larger range 
of Reynolds numbers, and for a variety of surface 
roughness patterns. So far, this demonstration is 
incomplete. Clearly, this tracking device must be built 
and tested. A definitive set of criteria for selecting the 
control law parameters, one that includes the effects of 
the probe hardware, must be developed. However, the 
results presented in this paper are very encouraging. The 
simulation results clearly confirm that the tracking 
probe can be designed to accurately measure the skin 
friction in a real-time environment. 


13 

American Institute of Aeronautics and Astr onautics 




Figure 1. Subsonic correlation of base and viscous forebody drag coefficients. 



Viscous forebody drag coefficient 
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Figure 2. Visualization of the drag bucket. 
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Figure 3. Conceptual schematic of boundary layer tracking probe. 
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O Closed-form solution, eqn. 12 
□ Edge solution, eqn. 8b 
O Nonlinear solution, (5 = 0 
A Nonlinear solution, (3 = - 0.25 Eqns. 4, 5, 
K Nonlinear solution, (3 = 0.25 ant * ® 

A Nonlinear solution, (3 = 1.0 


Note: Skin-friction coefficient 
referenced to wetted area 



Figure 5. Effect of pressure gradient on the law of the wake solution methods. 
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Figure 6. Effect of surface roughness on the law of the wake solution methods. 



Figure 7. Schematic showing location of boundary layer instrumentation on A-5A aircraft. 
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Figure 8. Schematic showing location of boundary layer instrumentation on XB-70-1 aircraft. 



Figure 9. Linear-log plot of boundary layer rake data. 
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Figure 10. Comparison of law of the wake model 
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Figure 1 1. Comparison of law of the wake model against XB-70-1 wing surface boundary layer rake data. 
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Figure 12. Summary of flight data, skin-friction coefficient comparisons. 
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Figure 13. Effect of the control gain parameter on probe response. 



Figure 14. Velocity profile exponent. 
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Figure 15. Response of probe control law to noisy boundary layer. 
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Appendix A: Integral 
Law-of-the-Wake Analysis 

This analysis develops an integral relationship that 
relates the longitudinally averaged skin-friction 
coefficient to a local measurement of the boundary layer 
thickness. The analysis is based on Coles’ empirical 
“law of the wake.” The law of the wake relates the outer 
(turbulent) boundary layer variables to the 
nondimensionalized boundary layer coordinates and can 
be written in terms of the velocity defect as 



90 

transitional. Spalding has offered a transitional 
formula for this buffer layer; but with minimal loss in 
accuracy a simpler model can be developed by equating 
A- 1 and A-2 


1 — / — V 

a/ J trans 


5 

4.878Ilcos2 

trans 

V 2 


_2 5 


(A-3) 


-2.439 In 


r>7 


Allowing that the cosine term in equation A-3 is 
approximately 1, equation A-3 may be solved 
numerically, to produce 


5. 

4.878TC 1-sin 2 

~ny~ 

V 2 

L 1 

_28_ 


2.439 hi 





4.8781Tcos 


~ny~ 

-2.439 In 

AT 

_28_ 


.81 


(A-l) 


In equation A-l, n is the wake parameter and is 
proportional to the local pressure gradient. Equation A- 1 
is a general expression that is valid outside of the 
laminar sublayer for plates, curved surfaces, surfaces 
with longitudinal pressure gradients, and rough 
surfaces. For a flat plate with zero pressure gradient, the 
wake parameter has an approximate numerical value of 
n ~ 0.55. 


y»ans~ 10 -81 (A-4) 

Thus for y + < 10.81 equation A-2 is used as the 
velocity defect model and for y >10.81 equation A-l 
is used as the velocity defect model. 


Momentum and Displacement Thickness Calculation 

The fundamental definitions for the (incompressible) 
momentum and displacement thickness are 



(A-5) 



(A-6) 


Correction for Laminar Sublayer 

Outside of the laminar sublayer, equation A-3 has 
been found to correlate well with measurements; 
however, within the laminar sublayer the correlation 
breaks down, hi the laminar sublayer the correlation is 
better represented by the relationship, u 1 = y + , and the 
resulting velocity defect law in the laminar sublayer is 


If the integrals of equations A-5 and A-6 are 
integrated across the depth of the boundary layer, using 
equation A-2 to evaluate the velocity defect within the 
laminar sublayer (y < 10.81 ) and equation A-l in the 
turbulent layer (y > 10.81 ), then an integral form of 
the law of the wake can be derived. The laminar 
contribution to the momentum and displacement 
thickness are given as 


1- 



= 1 - 





(A-2) 


Equation A-2 strictly holds for y ' < 5 and equation 
A- 1 strictly holds for y > 25 . Between these limiting 
values, the boundary layer velocity profile is 


0. = 58.428-- 

lam R e 


1 - 7.207 



(A-7) 


^ lam 


58.428 


Re,. 


(A-8) 


®turb 



|l.72463[II + 


4615411 + 9.45545 II + 5.9486] 


(A-9) 
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The turbulent contribution is given as 


Relationship of Wake Parameter to Local Pressure 
Gradient 



{ 1.72463 [IT + 1]} 


(A- 10) 


Adding the laminar contributions and turbulent 
contributions together the total normalized momentum 
and displacement thickness equations become 


l = (~ x »Ff){ L ? 2463[n + l]} 

_ C 7 [4.46154n 2 + 9.45545 n + 5.9486] 

V 1 (All) 

+ J-[ 32.0797- 26.3483 n] 

Re x 

+ 68.208,/^“ [n-0.978711][n + 3.09803] 

and 


As mentioned earlier in this appendix, the wake 
parameter (II ) is related to the local pressure gradient. 
The wake parameter has been correlated to the local 
pressure gradient by White 20 as a polynomial curve fit 
of the form 


dp 

* — 

0.42II 2 + 0.76n - 0.545 = R = — (A- 13) 

c f q 

In equation A- 13 (3 is the Clauser parameter. Solving 
for II hi terms of (3 and substituting into equations 
A- 11 and A- 12 the momentum thickness and 
displacement thickness equations can be written with 
the nondimensional pressure gradient as an explicit 
parameter in the equations 


8 


* 


X 


C?V^)(l-72463[n+ 1]} 
+ [32,0805 -26.34831T] 


(A- 12) 


0 

x 


[0. 16425 + 2.661 1670.88881 + (3 ] 
[-[10.4874+ 10.6227(3 + 2.1327670.88881 + J 


55,9194 - 40.656370.88881 + [3 

+ Re x 

JcJ [- 137.422 + 162.4P + 32.6051 70-88881 + (3| 


Re. 


(A- 14) 


and 


— = - /c7 [[0.16425 + 2.6611670.88881 + B]] 

x xV J 

. (A- 15) 

55,9194-40.656370.88881 + P } 


25 

American Institute of Aeronautics and Astr onautics 



(B-3) 


Appendix B: Unsteady Boundary Laver Model 

In a real-world turbulent boundary layer there are 
random elements whose characteristics can only be 
statistically predicted. Mathematically, this random 
element is treated by decomposing the boundary layer 
velocity components into time-mean and perturbation 
components. 

u[x, y,z,t] = U\x,y,z], + u'[x,y,z,t] 
v[x,y,z,t] = V\x,y,z] + v'[x, y, z, t\ (B-l) 
w[x, y, z, t ] = W[ x, y, 2 ] + w[x, y, z, t] 

The point-to-point variations of the perturbation 
components are inherently intermittent and can never be 
theoretically predicted. The law-of-the-wake analysis 
presented in the paper deals only with the time averaged 
velocity components of the boundary layer. However, 
the dynamic flow components will effect the 
performance of the boundary layer tracking probe and 
this flow intermittency must be modeled in the probe 
simulation. In this paper the random components within 
the boundary layer are modeled using data taken from 
Klebanoff. 19 The intermittency factor is the fraction of 
time that the flow is turbulent at a given perpendicular 
distance away from the wall. The intermittency factor is 
plotted as a function of the normalized perpendicular 
distance in figure B-l. The raw data from Klebanoff ’s 
report are presented along with a curve-fit model 


u[x, y, z, f] = U[x, y, 2 ] + |«'|yn[f] 
y\x, y,z,t.\ = |v'|yr|[f] 
w[x,y,z,t] = |vt/|ynl>] 

In equation B-3 t|(7) is a uniformly distributed 
random number, and the terms \u \ , |v| , and |w| are the 
expected values for the perturbation magnitudes. These 
parameters are also derived from curve fits of data 
presented by Klebanoff. These data are presented in 
Figure B-2. The unusual shape of the data curves was 
difficult to curve-fit. Instead these data were loaded into 
the simulation as tables and their values were computed 
as a function of y/8 using a log-linear interpolation. 

The resulting model produces a boundary layer with a 
ragged shape. In the outer regions of the boundary layer 
there are turbulent fluctuations which extend beyond the 
normally defined boundary layer edge. Conversely, 
there are intermittent regions inside of the normally 
defined boundary layer where the velocity is equal to the 
free-stream velocity. At each time point the dynamic 
pressure profile within the boundary layer is computed 
by taking the inner produce of the velocity components 
in equation B-3 and using this result to scale the 
undisturbed dynamic pressure profile. 


unsteady 


1 + 


2 Uu' + 1 ( 2 + v' 2 + w' 2 " 

u 2 


4 steady (B ' 4 ) 


y = 



-1.4 


(B-2) 


This intermittency function is used as a probability 
density function which describes the tune distribution of 
a particular velocity component. Assuming that the 
time-mean y and 2 velocity distributions are 
approximately zero, then the intermittent boundary 
layer model is written as 


In equation B-4 u\ v’, and w’ are the random velocity 
perturbations evaluated using the formulas of 
equation B-3, and U is the steady longitudinal velocity 
in the boundary layer at location y. The intermittent 
boundary layer thickness is computed as the value of y 
where 

5 -> Ju 2 + 2Uu' + u' 2 + v' 2 + w' 2 = U e (B-5) 


26 

American Institute of Aeronautics and Astronautics 




Figure B-2. Fluctuating velocity components in the turbulent boundary layer. 


27 

American Institute of Aeronautics and Astronautics 




References 

1 Wong, Thomas J., Charles A. Hermach, John O. 
Reller, Jr., and Bruce E. Tinling, Preliminary Studies of 
Manned Satellites — Wingless Configurations: The 

Lifting Body, NACA Conference on High-Speed 
Aerodynamics: A Compilation of the papers presented, 
NASA TMX-67369, 1958. 

2 Hoemer, Sighard F., Fluid-Dynamic Drag, Self- 
Published, Brick Town, New Jersey, 1965. 

a 

Saltzman, Edwin J., K. Charles Wang, and Kenneth 
W. Iliff, “Flight-Determined Subsonic Lift and Drag 
Characteristics of Seven Lifting-Body and Wing-Body 
Reentry Vehicle Configurations With Truncated Bases,” 
AIAA Paper 99-0383, January 1999. 

4 Whitmore, Stephen A., and Timothy R. Moes, “A 
Base Drag Reduction Experiment on the X-33 Linear 
Aerospike SR-71 Experiment (LASRE) Flight 
Program,” AIAA 99-0277, January 1999. 

5 Miller, Linda, Michael Hect, Martin Buehler, Amin 
Mottiwala, Beverly Eyre, and hidrani Chakraborty, 
“Ceramic Hybrid Electromechanical Systems 

(CHEMS),” NASA Tech Briefs, Vol. 23, No. 2, February 
1999, pp. 61-62. 

6 Coles, Donald, “The law of the wake hi the turbulent 
boundary layer,” Journal of Fluid Mechanics, Vol. 1, 
1956, pp. 191-226. 

7 Garringer, Darwin J., and Edwin J. Saltzman, Flight 
Demonstration of a Skin-Friction Gage to a Local Mach 
Number of 4.9, NASA TN D-3830, Feb. 1967. 

8 Quinn, Robert D., and Frank V. Olinger, Flight- 
Measured Heat Transfer and Skin Friction at a Mach 
Number of 5.25 and at Low Wall Temperatures, NASA 
TM X-1921, November 1969. 

9 Naughton, J. W., and J. L. Brown, “Surface 
Interferometric Skin-Friction Measurement Technique,” 
AIAA Paper 96-2 183, June 1996. 

10 Hopkins, Edward J., and Earl R. Keener, Study of 
Surface Pitots for Measuring Turbulent Skin Friction at 


Supersonic Mach Numbers — Adiabatic Wall, NASA 
TN-D-3478, July, 1966. 

1 1 Keener, Earl R., and Edward J. Hopkins, Accuracy 
of Pitot-Pressure Rakes for Turbulent Boundary-Layer 
Measurements in Supersonic Flow, NASA TN D-6229, 
March 1971. 

12 Mills, Anthony F., Heat and Mass Transfer, Irwin 
Publishing Co., Homewood, Illinois, 1992. 

13 Clauser, Francis H., “Turbulent Boundary Layers in 
Adverse Pressure Gradients,” Journal of Aeronautical 
Sciences, Vol. 21, No. 2, 1954. 

14 White, Frank M., Viscous Fluid Flow, McGraw- 
Hill, Inc., New York, 1991. 

15 Mills, A. F. and Xu Huang, “On the Skin Friction 
for a Fully Rough Flat Plate,” Journal of Fluids 
Engineering, Technical Briefs, Vol. 105, September 
1983, pp. 364-365. 

I6 Rade, Lennart, and Bertil Westergren, Beta 
Mathematics Handbook, 2nd ed., CRC Press Inc., Boca 
Raton, Florida, 1990. 

17 Fisher, David F., and Edwin J. Saltzman, Local Skin 
Friction Coefficients and Boundary-Layer Profiles 
Obtained in Flight from the XB-70-1 Airplane at Mach 
Numbers up to 2.5, NASA TN D-7220, June 1973. 

18 Saltzman, Edwin J., and David F. Fisher, Some 
Turbulent Boundary-Layer Measurements Obtained 
from the Forebody of an Airplane at Mach Numbers up 
to 1.72, NASA TN D-5838, June 1970. 

19 Klebanoff, P. S., Characteristics of Turbulence in a 
Boundary Layer with Zero Pressure Gradient, NACA 
Report 1247, 1955. 

20 Spalding, D. B., and S. W. Chi, “The Drag of a 
Compressible Turbulent Boundary Layer on a Smooth 
Flat Plate With and Without Heat Transfer,” Journal of 
Fluid Mechanics, Vol. 18, 1964, Cambridge University 
Press, New York, New York, pp. 1 17-143. 


28 

American Institute of Aeronautics and Astronautics 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, 
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, 
VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 

January 2000 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


4. TITLE AND SUBTITLE 


5. FUNDING NUMBERS 


A Real-Time Method for Estimating Viscous Forebody Drag Coefficients 


6. AUTHOR(S) 

Stephen A. Whitmore, Marco Hurtado, Jose Rivera, and Jonathan W. 
Naughton 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Dry den Flight Research Center 
P.O. Box 273 

Edwards, California 93523-0273 


WU 529-50-04-T2-RR-00-000 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


H-2394 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NAS A/TM-2000-2090 1 5 



11. SUPPLEMENTARY NOTES 

Presented as AIAA 2000-0781 at the 38th Aerospace Sciences Meeting and Exhibit, January 10-13, 2000, Reno, Nevada. 
Stephen A. Whitmore and Marco Hurtado, NASA Dryden Flight Research Center, Edwards, CA; Jose Rivera, University 
of South Florida, Tampa, FL; and Jonathan W. Naughton, University of Wyoming, Laramie, WY. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified — Unlimited 
Subject Category 05 

This report is available at http://www.dfrc.nasa.gov/DTRS/ 


13. ABSTRACT (Maximum 200 words) 

This paper develops a real-time method based on the law of the wake for estimating forebody skin-friction 
coefficients. The incompressible law-of-the-wake equations are numerically integrated across the boundary 
layer depth to develop an engineering model that relates longitudinally averaged skin-friction coefficients to 
local boundary layer thickness. Solutions applicable to smooth surfaces with pressure gradients and rough 
surfaces with negligible pressure gradients are presented. Model accuracy is evaluated by comparing model 
predictions with previously measured flight data. This integral law procedure is beneficial in that skin-friction 
coefficients can be indirectly evaluated in real-time using a single boundary layer height measurement. In this 
concept a reference pitot probe is inserted into the flow, well above the anticipated maximum thickness of the 
local boundary layer. Another probe is servomechanism-driven and floats within the boundary layer. A 
controller regulates the position of the floating probe. The measured servomechanism position of this second 
probe provides an indirect measurement of both local and longitudinally averaged skin friction. Simulation 
results showing the performance of the control law for a noisy boundary layer are then presented. 


14. SUBJECT TERMS 


Base drag, Boundary layer, Hypersonic vehicles, Skin friction 


17. SECURITY CLASSIFICATION 
OF REPORT 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 


Unclassified 


Unclassified 


Unclassified 


15. NUMBER OF PAGES 

34 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 

Unlimited 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 


NSN 7540-01-280-5500 




























